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Abstract 

We evaluate two-loop b ^ s'-f matrix elements of all the four-quark operators containing 
no derivatives. Contrary to previous calculations, no expansion in the mass ratio mc/mb is 
performed, and all the possible Dirac and flavor structures are included. Consequently, we are 
able to provide the last item in the NLO analysis of 5 — »• Xgj that has been missing so far, 
namely the two-loop matrix elements of the QCD-penguin operators. Due to smallness of the 
Wilson coefficients of those operators in the Standard Model, their effect on the branching ratio 
is small: a reduction by roughly 1%. We find BR[B Xs-fU^yis cev = (3.57 ± 0.30) x 10"^ 



1 Introduction 



The decay B — > Xs'j constitutes a stringent test of the Standard Model (SM) and many of its 
extensions. In the theoretical prediction for its branching ratio, crucial role is played by the 
NLO QCD corrections to the 6 — 57 partonic amplitude. In the introduction to our previous 
article |1|], status of the NLO calculations has been summarized. The only missing elements 
were two-loop matrix elements of the so-called QCD-penguin operators. 

The present paper is devoted to evaluation of these matrix elements. Our results for two- 
loop diagrams are presented in such a manner that they can be applied to an arbitrary extension 
of the SM where additional four-quark operators arise, e.g. to the generic MSSM. Thus, apart 
from completing the NLO QCD calculation in the SM, we provide an important ingredient for 
analyses of new physics effects. Moreover, contrary to previous calculations |0, |l|, no expansion 
in the mass ratio mc/mh is performed. 

The article is organized as follows. In the next section, the relevant definitions are collected. 
In section ^, we summarize our results for the matrix elements, and describe consequences for 
BR[i? — ^ Xg'-f] in the SM. Section |^ contains details of the calculation for one set of the two- 
loop diagrams. For the remaining sets, we present only the final expressions in section ^ The 
renormalization procedure is described in section ^. Section |^ contains our conclusions. 



2 The effective Hamiltonian 

In the SM, the 6 — > 57 transition is mediated by the effective Hamiltonian^ '0 

n,s = — y=-v*v,,Y.CkPk, (2.1) 

where Ck are the Wilson coefficients and Pk stand for the following operators: 
Pi = {sl1,T''cl){cl1^T%l), 

P2 = {sLli,CL){cLl^hL), 

P3 = {sLl^,bL)J2(^^^(l)^ {q = u,d,s,c,b) 

g 

P4 = {sLl,T%L)T.i<l^'T^<l)^ 
g 

P5 = {sLli,7ulpbL)Y.i^^''^''^''(l)^ 



^ It is written in terms of bare quantities here. The renormahzed interaction terms can be found in eq. (6.1). 
^For simphcity, we set Vub to zero in aU our formulae. However, its non-zero value will be included in the 
phenomenological analysis, as in eq. (3.7) of ref. H. 
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In order to match our computation with the NLO results for Ck, we adopt here the same 
operator basis as in ref. @]. 

Instead of the original Wilson coefficients Ck it is convenient to use certain linear combina- 
tions of them, the so-called "effective coefficients" 

Ck, for A; = 1, 6, 

C7 + ELi2/.a, for A; = 7, (2.3) 

C8 + E■=lZ^Ci, for A; = 8. 

The numbers i/i and Zi are defined so that the leading-order b —>■ s'j and b —>■ s gluon matrix 
elements of the effective Hamiltonian are proportional to the leading-order terms in Cf^ and 
C^ , respectively. In the NDR scheme, 
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The renormalization group equations for the MS-renormahzed effective coefficients read 

^^-^cfi^^) = cfif^h]fif^). 

Explicit results for the coefficients 

Cfif^) = ci'^^'i^^) + 
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can be found in ref. (see also eq. ( |6.3| ) here). The partonic decay rate from eqs. (29)-(32) 
of that article can be expressed as follows: 



where, for /i ~ nib, 
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Figure 1: Two-loop IPI contributions to the matrix elements of 



For k = 1, .., 6, the quantities in eq. (|0| ) can be found by calculating the two- loop 6 — »• S7 
matrix elements of the four-quark operators P^. The relevant Feynman diagrams are presented 
in fig. m, where the square boxes denote the operator insertions. On the other hand, rg is given 
by the one- loop matrix element of Pg. 



3 Final results for rk and consequences for BR[i? ^sl] 

Calculation of in eq. ( p.9|) has been the main goal of the present paper. For k = 1,2, 8, we 
confirm the findings of refs. 0, For k = 3, 6, our results are new. Altogether, they read 
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f|-|[a(^)+6(^)] + ||zvr. 
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+ 2[a{z) + h{z)]-f^m 



+ j^, + 'iX,-a{l) + 2h{l) + 

1 - ^ - i^^ + + W) + 2&(^) - W.^^ 



(3.1) 



56680 + |2. ^ v^x, - 16a(l) + 326(1) + ^^vr. 



27" 

243 ' ~'~ 9 - — V-/ I —V-/ ' 81 

W - g| - i^^ - f + f Kl) + 12a(z) + 206(z) - ^^vr, 

_ 44 _ 8 2 I 8 • 

rg — -g- — 277r + givr, 

where z = ml/ml. Similar ratios for the light quarks {u, d, s) have been set to zero. The 
constant is given by 

Xb= dx C dy dvxy\-a[v + x{l-x){l-v){l-v + vy)]c^-Q.lQU. (3.2) 
JO Jo Jo 
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Figure 2: Functions a(z) and b(z). Dotted lines represent their expansions at 2; = up to 0{z^). 
Exact expressions for the functions a{z) and h{z) are as follows: 



a{z) 



dx dy dv {[2 — V + xy{2v — 3)] ln[t>2; + x(l — x){l — v){l — v + vy)] 



+ [1 — V + xy(2v — 1)] \n[z — ie — x{l — x)yv]} + ^ + q^vt 



(3.3) 



h{z) 



ln2; + + 81 ^ 
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Vl^zf{z) - y [Iz - 1) f{z)' 



1 \y\y^ - l)x{l - x) + {2 - y)ui In Mi + {2y^ - 2y - l)u2 In U2 



(3.4) 



with Mfc = y^x{l — a;) + (1 — y)z, 



-1 = +i and 



/(^)=^^(l-4^) (ini^ 



/l=4l 



2i^(4z-l)arctan^7^. (3.5) 

Additive constants in the functions a and b have been chosen in such a manner that a(0) = 
6(0) = 0. Our results for a(l) and 6(1) read 



a(l) ~ 4.0859 + -i7r, 
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There is no need to apply numerical integration for z = m^/ml ~ 0.1, because, as illustrated 
in fig. ^ both functions are then accurately given by their expansions in z 
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Table 1: Real and imaginary parts of for two different values of mc/rrih- 

The numerical values of for two different values of rric/ mi, are presented in table |I|. We 
observe that the real parts of and are considerably larger than the remaining ones. If r4 
were as large as r^, its effect on the SM branching ratio would be around 6%. On the other 
hand, the size of is somewhat artificial — it is due in part to the lack of a "natural" factor 
^ in the definition of P5 (U)- 

In order to calculate BR[i? — > X^'y], we use the NLO formulae collected in ref. 0. The only 
elements that must be modified due to non-zero values of r^, rg are the so-called "magic 
numbers" . Their updated values are presented in table 0. Once they are used, and the same 
computer program with the same numerical inputs is applied,^ we find 

BR[B ^ Xs^Te^'^^TLv ^' = (3.57 ± 0.30) x 10-^ (3.10) 

which is only a little lower than (3.60 ± 0.30) x 10""^ in eq. (4.14) of ref. For the cutoff 
energy of ^rrib (which corresponds to the fake "total rate"), the effect on the branching ratio 

One exception is Sew from eq. (4.6) of ref. ^ which we change to 0.0071 according to ref. Q. This 
modification shghtly diminishes the net effect on the branching ratio. 
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Table 2: Update of table 1 from ref. 0. Only 4, dt, dl, dl"" for k = 3, ...,8 are affected. 



is the same: the central value decreases from 3.73 x 10^^ to 3.70 x lO"'^. 

Such a small effect of the non-zero values of r^, tq is due to smallness of the Wilson 
coefficients of the corresponding operators. In ref. 0], this effect was estimated to be around 
1%. For that reason, the NLO QCD computation was called "practically complete" already at 
that time. However, it becomes strictly complete only now, when we learn exactly what the 
effect is and in what direction it acts. 



4 Example of a two-loop diagram calculation: 

In the present section, we describe in detail the calculation of the two diagrams denoted by G3 
in fig. |1|. We choose {s l"^ ^CL){cL'~i^hL) for the inserted operator. Here, contrary to refs. [jl], 0, 
no expansion in z is performed, because the limit z ^ 1 is relevant for operators that con- 
tain three 6-quarks. The UV divergences are regularized dimensionally. We treat 75 as fully 
anticommuting in D = 4 — 2e dimensions. The s-quark is assumed to be massless. 
The considered Feynman integral can be written as follows 

G« = / fe2(f+fc)2 ^'7^(^ + ^)7pPl J,.eVPL6, (4.1) 
where 

{\p+}^+mc]'yu\^+mc\'yf,[p+/+mc] - \^+/-mc]'y^\f-mc]-fu\^+)^-mc]} (4.2) 

and 

A = [(j) + kf- ml] [{p + rf- ml] - m^] . (4.3) 
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Here, r and q stand for the outgoing photon and s-quark momenta, respectively. The momen- 
tum p runs inside the c-quark loop, while k is the gluon momentum. The polarization vector 
of the photon is denoted by e^. 

Let us first consider the one-loop subgraph J^^, for arbitrary off-shell momenta of the quarks, 
the gluon and the photon. The necessary one-loop integrals read 



J d^p/A = a, 
j d^p Pa/ A = bka + era, 
j d'^p PaPp/A = dmlgap + ekakp + fvaVp + g^kaTp + Takp) . 



(4.4) 
(4.5) 
(4.6) 



There is no need to consider PaVpP-y iii the numerator, because f^vf^p^f—f^^if^vf = p'^ipl^lu — 
Ivlnf}- The coefficients a, ...,g are not all independent. The following relations can be found 
by considering various contractions of the tensor integrals: 



b = 
d = 



la-c-2g-{c + 2f)^, 



(4.7) 



Semi 



a(A;2 - 4ml) + M^-r + k^ - 3r^) + 4f{k-r - 2r^) + 4g{k'^ - 2k-r)\ , (4.8) 
k-r + r"^ 



e = ^a + ^c + ^ + (c + 2/) 



(4.9) 



After performing the Dirac algebra with the above relations taken into account, one arrives at 
the following result: 



= + + fK + ^^/t + \aP;, + (c + 2f)P;i + gP^^^ 



(4.10) 
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(4.12) 



The structures are obtained from X^J!^ by interchanging /i ^ z/ and k 



r, I.e. 



Xj^\k,r) = Xj)^(r, k). Note that each of those structures vanishes under contractions with k'' 



and r^. This is a manifestation of the Ward identities k'^J, 



0. 



Once we have found the one-loop integral J^,^, we should complete the remaining Dirac 
algebra in G^^^ ([4.1|). Given our explicit expressions for the structures Xj^f^ and Xj^^, it is easy 
to verify that 

1pPlJ,.YPl = (2 + 2e)Pn \bV' + cK^ + fV^ + gV^, 



(4.13) 



Consequently, in our present calculation of the matrix element of {sL'~^^CL){cL'~^^hL), the terms 
proportional to in J^;^ ( [4.10| ) are irrelevant. Nevertheless, we have presented them explicitly 
here because they matter for other operators. 

From now on, we shall impose the on-shell conditions for the quarks and the photon: 
= = 0, q-r = = and (^ + f)h = m^jh. For an arbitrary scalar func- 

tion F{k'^, k-r, k-q), we find the following identities: 



j dPk F{k'^ , k-r, k-q)s'y'^{(^ + ^)V^^e^6 = [terms proportional to s^b and {e-r)sb] 



+mi,s(ff - 

j d^k F{e,k-r,k-q)sY{^ 

J d^k F{k\k-r,k-q)sY{^ 

J d^k F{k\k-r,k-q)sY{^ 
rribS 



{k + q) 



,k-r 



e 1 + 



2k-r 



mt 



J d^k F{k'^,k-r,k-q) 
+ ^)y,n.e'^b = [terms proportional to s^b and {e-r)sb], 



+ ]^)V^^e'^b = [terms proportional to sf^b and (e-r)s6]. 



(4.14) 

(4.15) 
(4.16) 



)V',,e^b 



d'^k F{k\k-r,k-q) 



= [terms proportional to s^b and {e-r)sb] 
2e(. + ,)^^-.^(l + (2-3e + 2e^)|;) 



(4.17) 
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Consequently, G^^^ ( |4.1| ) takes the following form 

G.f ^ = AT{l + 2e)TT^m^^' sPr {R^^^rribW -/^) + R^^^^rrib e-r + i?3^V } b, 
where 



(4.18) 



(1) 



g{k'^, k-r 



1-e 



6(A;^, A;-r) 



(A; + g) 



2^-^ ,2 

mi 



2k-r' 



2e(fc + g)'^ - A;M 1 + (2-3e + 2e^ 



mt 



As discussed in the next section, we do not need to calculate R^'^^ and R^,^'^' ■ The overall 
normalization factor in eq. (^4.18 ) is chosen in such a manner that the 0{1) part of R^^^ is equal 
to the quantity M3 from our previous article (see eqs. (2.6) and (2.11) there). 

The terms proportional to {k + g)^fc ■ r in the curly bracket of eq. ([4.19|) give vanishing 
integrals over k. Consequently, R^^^ simplifies to 

d^k 



mf 



2k-r' 



mt 



'(1) 



(4.19) 



4^^ 



:i + e)mf 



2r(l + 2e)7r^y {k + q)^ + ie 



1 + 



2k-r 



mt 



1 + (2-3e + 2e^) 



2k-r 



mt 



1 



(4.20) 



Now, we need explicit expressions for b and g. They can be easily found from eqs. ( [4.5|) and 
with the help of Feynman parameters 



6(A;^, k-r) 



g{k'^,k-r) 



m^'^Til + e) / dx I dy 

'Jo Jo ^ \-(k- 



X "(1 — x)~ 



ITX 



^/2r(l + e) / dx dy 



1 



X 



X] 



xy). 



lo Jo [_(fc_^^)2 + £!^]l+e 

Introducing another Feynman parameter and integrating over k, one obtains 



(4.21) 
(4.22) 



R 



(1) 



l + e 



4e(l - 



' dx dy 
Jo Jo 



dv 



(1 - e)v + xy[{l -v){2-3e + 2e^) - 1] 



(4.23) 



X "^(1 — x) v*^ [z — ie — x(l — x)y{l — v)]' 

where z = m'^/mf. After expanding the integrand to 0{e) and performing the easy integrations, 
we find 

(U 111 /"i /"i /"i 

Rl' = \ h-/ dx dy dv[l-v + xy{2v-l)]\n[z-ie-x{l-x)yv]. (4.24) 

8e 8 2 Jo Jo Jo 

In the next section, analogous results for all the diagrams in fig. |1] will be presented. 
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5 Results for the unrenormalized two-loop diagrams 



In the present section, we shall be interested in the unrenormalized two-loop IPI on-shell 
b s'j diagrams with insertions of the following 4-quark operators: 

1 

Q^"^ = -^(sPR7ai-7.„c)(c7""...7"^PK&). (5.2) 



Q^"^ = — T(^Pfi7ai-7«„c)(c7""...7"^PL&), (5.1) 



The corresponding diagrams in fig. |I| will be denoted by G^"'' and G^f^\ respectively. Their 
overall normalization is assumed to be the same as in eq. ([4 .11), e.g. 



= ^ / p(f+\)2 ^'7-(^ + l^)PRla. . . . laj,.e^l-- . . . (5.3) 

The operators Q^^\ Q^^\ Q^^\ Q^^\ Q^^) together with their mirror copies and analogous 
operators with color-octet currents form a complete set of dimension-six (sc)(c6) operators in 
four spacetime dimensions. Extending our results to Ai? = —AS* = 1 four-quark operators 
with other fiavor contents is straightforward, as we perform no expansion in the ratio nic/mi,. 

For the operators Q only the diagrams G'f'\ C'f'^ matter. The remaining ones in fig. |l] 
vanish due to chirality conservation in the charm-quark loop and QED gauge invariance. 

For Q("), the set 

gives no contribution on-shell, i.e. it cancels with the correspond- 
ing IPR diagrams, because it is proportional to the matrix element of a two-quark operator 
sPuh that vanishes by the equations of motion, up to a total derivative 



sPnh = -— 



(5.4) 



The diagram Gg"^ is just a product of two one-loop diagrams. It contains an IR divergence 
that cancels out only after adding the corresponding bremsstrahlung corrections. However, for 
all the Q*-"^ except Q^^^ (which is of no interest in the SM), the effect of those diagrams can 
be taken into account by the standard replacement — ^ Cf^ (see section P). Thus, we shall 
ignore Gg"^ in this section, i.e. we shall consider only G^""*, G^^^ for the operators Q*-"'^ 

By analogy to eq. ([4.18|) , we write 

g1") = 4r(l + 2e)7r^m,-4^ sPr {Rt^m,{/f - //) + R^'^^m^ e-r + i?^("V } ^ (5.5) 

and similarly for G^"^ Below, we shall present our results for the coefficients in front of (//—//) 
only. Due to QED gauge invariance, the remaining Dirac structures must cancel out when all the 
IPI and IPR diagrams are included. No explicit calculation of the IPR diagrams is necessary, 
as they can give no — //) structure. 
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We have found the following expressions for rP and 
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I I fl rl rl 
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+ 3^ In ^ _ n±6 _ 2 /^i rf^ j^i rf^ /^i dv xy In , ^ odd. 



(5.10) 
,(5.11) 

(5.12) 



n—2 I 3— n 
4e ~'~ 2 



Y(-^ + J,' dxJ.Uy J,' dv{'-^ In 



zv+x{l—x){l—v){l—v+vy) 



^Inz}, n 



(5.13) 



even, 



1 5 /3 TT^ ^ 

2 In 2 he h51n2 + 21n2z 

e 2 U 6 



(5.14) 



where F^") = + ] {\nz + \) - In'' z + ^ - '-^^^ and Wk = ^{ffg^- The variables 
and the function f{z) are as in b{z) in eq. ( p.4|) . 

We have presented R^^^ and -R4"'' for arbitrary n, while n = 1 was assumed otherwise. In 
those cases, results for other n can be found from the following simple relations: 



G 



(n+l) 



G 



(n+l) 



-1+e r^{n) 
2 , 


for k = 


1,2, 


(because 'jalij.l" = ( 


-2 + 2e)7^ 


+ l+e 

2 "^fc , 


for k = 


3,4, 


(see eq. (4.13)), 




2 '-^fc ' 


for A; = 


1,2,5, 


(because '~iaO'^v'^°' = 


-2ea^y). 



(5.15) 
(5.16) 
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The overall factor of e in the r.h.s. of the latter equation is the reason why we have calculated 
the C(e) parts of r'^\ R^^^ and R^^K They are necessary for the 0{1) matrix elements of Q^^\ 
i.e. they may matter beyond SM. 

For the color-octet analogues of Q^'^) and the expressions for i?^"^ and i?^"'' get multi- 
plied by additional color factors of — I (for A; < 4) and | (otherwise). Of course, the values of 
and ^i"^ remain the same for the mirror copies of Q^"-* and Q^"\ provided Pr is replaced 
by Pl in eq. (|5.5| ). 

Thus, apart from the trivial case of the present section contains complete results for 
the unrenormalized two-loop b ^ s'y diagrams with insertions of four-quark operators, in the 
SM and beyond.]] 



6 The MS-renormalized amplitude in the SM 



In the present section, we shall restrict ourselves to the SM operators Pk ( p.2|) and use eqs. (|5.6|) - 
( |5.14| ) to evaluate the MS-renormalized on-shell b s'j amplitude at the NLO in QCD. From 
this amplitude one reads out the coefficients r^, i.e. our main results presented in section ||. 
The renormalized effective Lagrangian readsQ 



eff 



c 



QCDxQED 



V2 



H \ C7{fi)Zf7Z^ZmP7 + Cs{fi)Z^Zm ( ZsgZgZi^Pg + ZsrP-, 



k=l 



The evanescent operators E^^^ and E^''' are as in the appendix of ref. ^ 



6 4 

ZkjPj + X] ^kjEj ) + ZkjZ^ZmPj 

(1) 



+ ... 



(6.1) 



e: 



(1) 



(1) 



(^L7m7M2 7M3 7M47M5&L) E,(g7'^^7^^7'^^7'^n^^g) - 2OP5 + 64P3, 
(^l7mi7m2 7m3 7m47m5^"&l) j:,{qr'r'r'r'r'T^q) - SOPg + 64P4. 



(6.2) 



The dots in eq. (|6.1|) stand for other evanescent operators that do not affect the NLO matrix 
elements. They are relevant for the NLO anomalous dimension matrix though. 

Below, in the calculation of the MS-renormalized amplitude, we shall use the MS-scheme 
renormalization constants, and implicitly assume that all the unrenormalized Z-loop matrix 
elements are multiplied by (47r)~''^e'^'''. The light {u,d,s) quark masses will be neglected. 

In the NLO calculation of 6 57 matrix elements, the renormalization constants Zg and 
Zg can be set to unity. As far as Z^ and Z^ are concerned, we need the one-loop expressions 



1 - and Z 

7re 



1 - 



Sire ' 



Up to operators that vanish in 4 spacetime dimensions, i.e. evanescent operators — see section ^. 
^ Note that the interaction terms in eq. (^.iD have opposite sign. 
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The renormalization constants Z^^ and Z^^ can be found from eqs. (20)0 and (45) of 
ref. 0. The only non- vanishing contributions to those constants at order 0{as) are 



A-Ke ' 



^63 - ana - 



The renormahzation constants can be found from the anomalous dimension matrix 7^'')'^*^ 



given in eq. (8) of ref. @][] 



7' 



7 



(O)eff 



-4 


8 
3 





2 

9 








208 

243 


173 

162 


12 








4 
3 








416 
81 


70 
27 











52 
3 





2 


176 
81 


14 
27 








40 
9 


100 
9 


4 
9 


5 

6 


152 
243 


587 
162 











256 
3 





20 


6272 
81 


6596 
27 








256 
9 


56 
9 


40 
9 


2 

3 


4624 
243 


4772 
81 




















32 
3 























32 
9 


28 
3 



The relation 
case, 



kj 



Skj + -^Ikj"^ holds for all the Z^j except (Zj^r, Zj, 



(6.3) 



'fc=l,...,6- 



In the latter 



a. 



IGvre 



(0) 

Tki 



J=3 



(6.4) 



7. 



eff 



.,6, 
.,6, 



(6.5) 



where 7 is related to 7''''^ by |^ 

7i7 + ELi Vkljk - Vjlrr - Zj'jsr, 

7j8 + EE=i Zk-fjk - Zj-f S8, 
7ji) 

with Ui and Zi given in eq. ( |2.4| ). The numbers yk and yj' parametrize the unrenormalized 
on-shell one-loop b ^ s'y matrix elements of the 4-quark operators 

2e 



when i = 7 and j = 
when i = 8 and j = 
otherwise, 



(-Pfc)l loop 



— ) ^yk{P7)tree + Oie'), 

rribJ 
rrib/ 



{Ej )lloop — 

An easy one-loop calculation gives yi 
y3 = Qd = -|, 
^5 = 4(5 - 3e)Qd, 



{yk 



lim?/fc) 



limyf) 



(6.6) 
(6.7) 



1/2 



16(4 - 25e)Qd, 



and 

y& 



4 
"9' 



QdCp - 

4(5 - 3e)g,Cp, 



(6.^ 



yX = 16(4 - 2ht)QdCF. 



^ In that equation, the expression "(counterterms due to ZV)'^ has been missed. 



Signs of 7^7^^^^ and 7^3^°" are fixed by our sign convention inside D^ijj = (d^ + igG°T" + ieQA^) tp. 



(O)cff 
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The unrenormalized one-loop matrix elements of P7 and can be expressed as follows: 

(i^7)lloop = ^f7(P7)tree, (6.9) 
(P8)lloop = (£^)''(-Z87 + ^r8)(P7)tree + 0(«',e). (6.10) 

No explicit expression for is necessary for our purpose. It is enough to mention that r-j is 

UV-finite (note that Zi-jZ,^Zm = 1 + 0(ag)), but it contains an IR divergence. We assume that 

this IR divergence is regulated with a small gluon mass. 

The quantity rg reads 

44 8 2 8 . ,^ ^ ^ , 

rs = TT H — m. (6.11) 

^ 9 27 9 ^ ' 

We have verified this result during the calculation of i^^^ and I^^^ . 

In the following, we shall need two unrenormalized matrix elements of evanescent operators 

(Pf+iPf-iQ« + iQ(°))2 loop = ^e4(P7)tree + 0(e), (6.12) 
(p--+ip--_8Q«+g(0))2i„,p = ^e6(P7)tree + 0(e), (6.13) 

where are just the cc-parts of the respective P^. Note that {P^'^)2 loop = {Pi^)2 loop = due 

to color factors. We find 64 = Qu — ^Qd and eg = IQQu — ^Qd- These results are independent 

of m^. Thus, after appropriate replacement of quark charges, they hold for any flavor circulating 

in the loop. Consequently, we can use eqs. (|6.12|) and ( |6.13|) to express all the two-loop diagrams 

with Dirac traces in (Pa;)2 1oop as linear combinations of the results from section ^, up to the 

following additive constants obtained from 64 and eg: 

8 32 
Da = uQu + dQd- —{u + d)Qd, and = 16{uQu + dQd) - —{u + d)Qd, (6.14) 

where u = 2 and d = 3 are the numbers of active up and down flavors, respectively. 

Such an approach is much more convenient than calculating diagrams with Dirac traces 
separately, because 64 and eg can be (and are) found using expansion in external momenta and 
computer algebra, as in the so-called matching computations (see e.g. ref. 0). More precisely, 
such an expansion is applied to diagrams with subtracted sub divergences, because subtracted 
matrix elements of evanescent operators must be local (i.e. polynomial in external momenta). 
Once they are found, we add the one-loop sub divergences again, setting the external momenta 
on shell. Then, the only relevant subdivergence is proportional to (P4)iioop that we already 
know. 

Now, we are ready to calculate the unrenormalized two-loop matrix elements of P^. Let us 
write them as follows: 

(Pa:)21oop = ^ (—X\xk + r-ryk){P7)tree + Oie), (6.15) 
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where the terms proportional to originate from Gq"""* diagrams. Evaluation of Xk amounts to 
forming appropriate linear combinations of eqs. ( |5.ti| )-( ]^7T^ ) with z ^ 0, z ot 1, according to 
the Dirac and flavor structure of Pk- Explicitly 

1 r, 
"6 



Xi 



Q,H^^{z) + Q^H\ 



34 



X2 

a;4 



QdH^^{z) + Q^Hl'^{z 



15 



- (i/ff (0) + H['^{1) + H['^) + CfQ.K 



r(l) 



(1) 



r(l) 
'55 



X5 



Q, [i/;?(o)+i/;?(i)+i/(') 
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+ Qu [HIM + ^34(^)] + Q<1 W[2{^) + 2i^34(0) + ii'[2{^) + ^m(I) 



(i/if (0) + Ef}{\) + + c^g.if. 



(3) 
55 



(6.16) 



where 



KM 



fc=i 



(n) 
A: 5 



(0), 



+ (20 - 12e)i?, 



(1) 



n 



0,1, 
3. 



(6.17) 



The relations between S'j"'' and (-R^"'') are the same as between ^^"^l) and (-Rj-"^) ^- The 

necessity of introducing S^"''' and 5*^"^ follows from the fact that ordering of the Dirac matrices 
in P5 and Pq is opposite to the one in Q^^^ and Q^^\ 

After all the above substitutions, one finds (up to o{e)) 



Xi 

X2 
X3 
X4 

Xq 



46 I 833 
243e ' 729 



f[a[z 



248 I 2932 8 „2 , Stt , 32 r, ( T \ J- Oh( T \ J- ^'^^ ott 



46 
"243e 



1031 I A„2 
729 ~'~ 81 



__l|X, + ia(l) + |6(l) + 26(^)- 



5552 _|_ 42424 
81e ~'~ 243 



4lT 

327r 



243 ' 



- + IS + ^^^^ - 16a(l) + 326(1) + 



2336 
81 



ITT, 



4588 
" 243e 



14378 
729 



80^2 1671 64 y 



9V3 27^ 



fa(l) + f6(l) + 12a(z)+ 206(2) 
15 



(6.18) 



3016 
243 



where the constant as well as the functions a{z) and b{z) have been already defined in 
section |[ 

Once we have found the unrenormalized matrix elements of Pi, Pg^ it is straightforward 
to calculate the renormalized ones. The MS-renormalized amplitude for b s'j decay reac 



M 



tree + (-^7)1 loop] + [{Ps)! loop + Zsj{Pj) 



tree I 



+ E + f^C'fV)! E [ZkjZ^ + (1 - 2e)(Z„ - 1)4,] (P,)i loop 

k=l L J I i=l 



+ (Pa:)2 1oop + Zk7{P7)tree + E ^fi (^^j'^ ) 1 loop ^ + 0{al € 



(6.19) 



i=3 



When comparing the above equation with eq. ( |6.1| ), one should remember about the identity 

The appearance of Z^p and Zm in eq. ( |6.19| ) can be understood without calculating any 
diagram. It is enough to remember that insertions of the Z^-counterterms on internal quark 
propagators always cancel with the Z^-counterterms at the ends of those propagators. Thus, 
we are left with a single power of Z^ that corresponds to two external quark lines. The term 
proportional to {Z^ — 1) is found from 



[(-Pfc)l loop] 



{Pk) 



1 loop; 



(6.20) 



where (Pfc)iioop is given in eq. (|6.6| ). Remember that (P7)tree depends linearly on mb(/x). 

When all the matrix elements and renormalization constants are substituted into eq. ( |6.19| ), 
one finds that all the 1/e singularities cancel as they should. Furthermore, when the effective 
Wilson coefficients are used, all the logarithms of m,f,/n are found to be multiplied by numbers 
from the 7th column of ■y^^^^^ (|6.3| ). Thus, the amplitude M takes the form 



1 , 



'^(O)efr tts ^(l)eff' 
^7 + -^^7 ^ 



I ^(O)efr 
+ ^ 1^ 

fc#7 



, (O)efr , 



(6.21) 



From the above expression, we have read out the results for Vk that have been already given in 
eq- (O- 

^ Up to an overall normalization factor of —^^V*gVtb- 
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7 Conclusions 



We have computed the two-loop b ^ s'y matrix elements of all the four quark operators 
containing no derivatives. This allowed us to complete the NLO QCD calculation of i? ^ 
decay by including for the first time the two-loop matrix elements of the QCD-penguin operators 
P3, ...,^6- The values of the corresponding parameters that enter the branching ratio are 
collected in table ^ They are also valid in extensions of the SM. 

The Wilson coefficients of QCD-penguin operators are small in the SM. Consequently, two- 
loop matrix elements of these operators affect the B X^'y branching ratio by around 1% 
only. In extensions of the SM, the Wilson coefficients might be larger, which would enhance 
their phenomenological significance. However, it should be emphasized that the modification of 
Ci with i = 3, 6 would not only have a NLO effect (via the two-loop contributions evaluated 
here) but also a LO one (via the effective coefficients in eq. ( |2.3| )). 

Since in certain extensions of the SM, new operators with different Dirac and color struc- 
tures may contribute, we have performed calculations that allow inclusion of such operators if 
necessary. The relevant results can be found in section |[ However, a complete NLO analy- 
sis in the presence of new operators would require extending the 3-loop anomalous dimension 
computation of ref. |Q, which is beyond the scope of the present work. 

On the technical side, occurrence of internal b-quark propagators required going beyond the 
techniques developed in refs. ^ that used expansions in z = ml/ml. Our exact calculations 
valid for any z show that the expansions performed in refs. |^ break down around z ~ 0.5, 
and are very accurate up to 2; ~ 0.3, i.e. well above the physical value of z ~ 0.1. 
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